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PROBABILISTIC VALIDATION OF HOMOLOGY COMPUTATIONS 

FOR NODAL DOMAINS 

By Konstantin Mischaikow-'^ and Thomas Wanner^ 

Rutgers University and George Mason University 

Homology has long been accepted as an important computable 
tool for quantifying complex structures. In many applications, these 
structures arise as nodal domains of real-valued functions and are 
therefore amenable only to a numerical study based on suitable dis- 
cretizations. Such an approach immediately raises the question of how 
accurate the resulting homology computations are. In this paper, we 
present a probabilistic approach to quantifying the validity of ho- 
mology computations for nodal domains of random fields in one and 
two space dimensions, which furnishes explicit probabilistic a priori 
bounds for the suitability of certain discretization sizes. We illustrate 
our results for the special cases of random periodic fields and random 
trigonometric polynomials. 

1. Introduction. The practical need to extract low-dimensional nonlin- 
ear structures from high-dimensional data sets has led to the introduction 
of topological methods in statistical analysis. There is a growing body of lit- 
erature [10, 13, 14, 26, 28] that addresses the problem of providing efficient 
algorithms for estimating topological properties of a fixed, but unknown, 
manifold X given a point-cloud data set that lies on or near X. The equally 
important complementary question is that of the accuracy of these estima- 
tions. To the best of our knowledge, the first result is due to Niyogi, Smale 
and Weinberger [25]. In this paper, the authors propose a stochastic algo- 
rithm for computing the homology of a given manifold X CM.'^ hy randomly 
sampling M points from the manifold and derive explicit bounds on the 
probability that their algorithm computes the correct homology. The prob- 
ability bound depends on the number M and a condition number 1/r. The 
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latter parameter encodes both local curvature information of the manifold 
X and global separation properties. More precisely, the inverse condition 
number r is the largest number such that the open normal bundle about 
X C of radius r is embedded in R*^ for all r <t. 

Although stimulated by the aforementioned results, the motivation for 
the work presented in this paper comes from the study of patterns produced 
by nonlinear evolutionary systems such as phase separation in materials, 
fluid flow and population biology. While the mathematical models for such 
systems are typically infinite-dimensional, understanding the time evolution 
of the multitude of often complicated patterns produced on an underlying 
two- or three-dimensional domain is of considerable, if not primary, interest. 
As is described in greater detail below, algebraic topology and, in particular, 
computational homology can be used to obtain new statistical measurements 
for either numerical or experimental data associated with these spatially and 
temporally complex systems [18, 19, 23]. As above, the accuracy of these ho- 
mology computations is an important question. However, because the data 
is obtained through numerical methods or by standard imaging techniques, 
the data points are not naturally modeled via a random sampling, but rather 
appear on a fixed grid. On the other hand, the patterns are constantly chang- 
ing and thus one is interested in the topology of a multitude of manifolds. 
Furthermore, at those points in time for which the topology of the pattern 
changes, the condition number 1/r becomes unbounded. 

These observations suggest the need for a complementary result concern- 
ing the accuracy of homology computations. To be more precise, consider 
a compact rectangular domain G := [a,b]'^ C M."^ and a smooth function 
u E C^(G,M). The nodal domains of u are given by the sub- and super-level 
sets 

iV± ■={xgG:±u{x)>0}, 

which can be interpreted as patterns on the domain G. These are the objects 
which we wish to identify through their homology, that is, we are interested 
in the statistics of H^[N^) as the functions u vary. 

Given the limitations of numerical and imaging (think digital camera) 
techniques, the nodal domains must be approximated. The following pro- 
vides a reasonable model for such an approximation. 

Definition 1.1 {Cubical approximation of a nodal domain). Consider 
a compact rectangular domain G := [a,b]'^ CM.'^ and a continuous function 
« : G — > M. Let M be an arbitrary positive integer. Define the equidistant 
M -discretization of [a, b] as the collection of grid points 

Xk '■= a + k ■ , for A; = 0, . . . , M. 
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Fig. 1. Nodal domain approximation using cubical sets in one dimension. 

The cubical approximations Q^j of the nodal domains = {x £ G: ±u{x) > 
0} are defined as the sets 

Qm := u| n [^^' + 1] • ' • • • ' ^d,k,) >0,ke{0,..., mY 

where k = (ki, . . . ,kd) and where Xjfl,...,Xj^M denotes the equidistant 
M-discretization of the jth component interval of G. For the case of one- 
dimensional domains, this definition is illustrated in Figure 1; for the two- 
dimensional case, see Figure 2. 
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At first glance, one might expect that for sufficiently smooth functions u 
and sufficiently large discretization size M, the homologies of the nodal 
domains and of their cubical approximations Q^j would coincide. Simple 
counterexamples to this exist; however, this assertion can be made correct 
from a probabilistic point of view. 

Consider a probability space {^},J-,¥) and let u:G x ^ M denote a 
random field over {Q,J^,¥) such that for P-almost all a; G fi, the func- 
tion u{-,io) :G — >M is twice continuously differentiable. We will impose the 
following two assumptions throughout this paper: 

(Al) for every x £ G, we have ¥{u{x) = 0} = 0; 

(A2) the random field is such that P{0 is a critical value of u} = 0. 

To simplify the notation, we follow the customary procedure of dropping 
the explicit w-dependence in the formulation of events, that is, we write, for 
example, F{u{x) = 0} instead of F{{u! G r2:?x(x,a;) = 0}). 

Theorem 1.2. Consider a probability space {Q,J-,¥), a domain G := 
[a, 6]'^ C and a random field x 17 — >R over P) such that for 

F-almost all uj G Q., the function u{-,u>) : G ^ M is twice continuously differ- 
entiable. For each lo £Q, denote the nodal domains of u{-,uj) by N'^{u)cG 
and given a positive integer M , denote their cubical approximations by Q^j{uj) 

If both (Al) and (A2) are satisfied, then for T-almost all u; G $7, the fol- 
lowing holds. For all sufficiently large values of M , the homology of the 
cubical approximations Q^,^{uj) matches that of the nodal domains N^{uj). 
In other words, there exists a random variable M : — > N such that 

¥{For all M>M one has if*(iV^) = H^Q^j)} = 1. 

Note that the random variable M is in general neither constant nor bounded. 

Remark 1.3. One can show that in the situation of Theorem 1.2, there 
actually exists a homeomorphism /i : G — > G which is isotopic to the identity 
such that h{N^) = Q^.j- In other words, the nodal sets and their cubical 
approximations are homeomorphic. Similar extensions are also valid for our 
probabilistic main results. Nevertheless, in view of the current state of ap- 
plications, we present all of our results within the homology setting. 

We do not provide an explicit proof of this result since it follows from the 
realization that for P-almost all w G 17, the level set w) = is a manifold 
with finite condition number 1/r < oo and a straightforward modification 
of the proof of [25], Proposition 3.1, then readily furnishes the above result. 

Observe that the hypotheses for Theorem 1.2 provide a fairly general 
framework in which to study the topology of random manifolds. However, 
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Fig. 3. Complicated microstructures generated from simulations of the Cahn-Hilliard — 
Cook model in two and three space dimensions. 

the motivation for this work leads us in a shghtly different direction. As an 
example, consider the problem of phase separation in binary metal alloys. 
During the last few decades, many models have been proposed for describing 
this process. Starting from the celebrated Cahn-Hilliard partial differential 
equation model [8, 9], numerous extensions and refinements have been pro- 
posed over the years. While all of these models are infinite-dimensional evo- 
lution equations they can be deterministic or stochastic partial differential 
equations and can contain nonlocal or other additional terms. Moreover, all 
of these models are phenomenological models in the sense that while they 
are based on important physical insights, they are not rigorously derived 
from first principles and so assessing the quality of a specific model becomes 
a central and crucial problem. Since the material microstructures which are 
generated by the phase separation process have direct implications for ma- 
terial properties, any such assessment should consider to what extent the 
complicated patterns by the models match experimental observations. 

In the above-described setting of evolution equations, patterns are usually 
described in a very specific framework. Returning to the materials example, 
the various models of Cahn-Hilliard type describe the temporal evolution of 
a phase variable u = u{t,x), where x £G describes the underlying container 
for the alloy. The scalar phase variable u describes the local alloy concentra- 
tion in the following way. If u is close to +1, the local concentration consists 
only of the first alloy component, values close to —1 represent the second 
component and values in between indicate corresponding mixtures of both. 
Typical patterns that are produced by the Cahn-Hilliard models are shown 
in Figure 3 and mathematical results have recently been obtained which 
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provide some insight into their dynamics [6, 7, 30]. Since we are interested 
in regions where one of the two components dominates, the microstructures 
are generaUy thought of as nodal domains of the function u, 

N^{t) = {xeG:±u{t,x)>0}. 

How weh do these nodal domains match actual experimental results? In a 
series of papers [20, 21, 24], Hyde et al. described experimentally obtained 
microstructures in iron-chromium alloys using a variety of techniques. One 
of their measurements was the handle number of the structure which, from a 
mathematical point of view corresponds to studying the first Betti number 
of the nodal domains N^{t) as a function of t. In [19], we used methods 
from computational homology to numerically study the Betti numbers of 
the nodal domains, in an attempt to unveil differences between two different 
models, the deterministic Cahn-Hilliard equation and the stochastic Cahn- 
Hilliard-Cook model. Since the underlying nonlinear evolutions were solved 
numerically, we do not have direct access to the nodal domains. However, for 
a fixed grid size bounded from below by the resolution used in the numerical 
approximation, the associated cubical approximation is easily obtained. 

This last comment indicates the limitations of Theorem 1.2. Ideally, given 
an evolution equation and a fixed cubical approximation of the nodal do- 
main, we would like to be able to estimate the amount of time for which 
the homology computations of the nodal domains are correct. This problem 
appears to be extremely difficult. Thus we pose the following, simpler, goal. 
Consider a random generalized Fourier series of the form 

oo 

(1) u{x,uj) = '^ak- gkii^) ■ '^k{x), u:G xQ^R, 

where G CW^ denotes a compact rectangular set and ip^, k gNq, denotes 
a sequence of basis functions which will be specified in more detail later. 
Assume that the gk are independent identically distributed real-valued ran- 
dom variables over a common probability space {^},J-,¥). The coefficients 
denote deterministic real numbers which will be related to the smoothness 
properties of the realizations u{-,uj):G ^M. Let M be an arbitrary posi- 
tive integer. Approximate the random nodal domains N^{uj) of the realiza- 
tion n(-,u;) by random cubical sets Q^j{io), as in Definition 1.1. Estimating 
the accuracy of our homology computations then reduces to the following 
problem: 

(P) Find sharp lower bounds on the probability 

F{H.iN^) = H,{Qt,)}- 

While this setting might seem restrictive, it does arise naturally in the 
context of evolution equations. If we consider linear stochastic partial dif- 
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ferential equations or linear deterministic partial differential equations orig- 
inating at randomly chosen initial conditions, then the solution at any later 
point in time will be a random Fourier series, as in (1). 

We would also like to point out that our setting is a special case of ran- 
dom fields. The geometry of random fields has been studied by a variety of 
authors; see, for example, [1, 2] and the references therein. These studies, 
however, almost exclusively consider the study of the Euler characteristic 
of nodal domains, which are called excursion sets in their setting. As was 
pointed out in [19], however, in many cases studying the Euler characteris- 
tic alone does not provide enough information. We therefore hope that our 
study of the Betti numbers in this paper will open the door to a more refined 
study of geometric properties in a stochastic setting. 

With these introductory comments in mind, let us now turn to the results 
of this paper. To focus on obtaining sharp lower bounds asked for in (P), we 
restrict our attention to one- and two-dimensional domains. 

Theorem 1.4. Consider a probability space P), a compact inter- 

val G = [a, b] in M and a random field ?/ : G x Q — > M over {^,J-, P) for which 
assumptions (Al) and (A2) hold and such that for ^-almost all G il, the 
function u(-,ti;):G^M is twice continuously differentiable. In addition, for 
X & G and 6 >0 with x + 6 £ G, let 

Paix, 6) = p|(7 • u{x,uj) >0,a-u^x + ^,uj^ < 0, 0" • u{x + S,U)) > 0^ 

and assume that there exists a constant Co > such that 

(2) Pa{x,6) <Co- 6^ for all a e {±1} and X eG with x + 6 gG. 

Then for every discretization size M , the probability that the homologies of 
N^{lu) and Q^_^{uj) coincide satisfies 

8Coib-af 



(3) F{H,{N^)^H,{Qt^)}>l 



3M2 



The proof of this theorem — which is fairly straightforward when viewed 
from the proper perspective — and its implications are presented in Section 2. 
Observe that the assumption in (2) lies at the heart of Theorem 1.4 and 
establishing its validity generally requires some work. In particular, deter- 
mining or estimating the constant Co has to be done with significant care 
since it has direct implications for the tightness of the resulting bounds. 
We obtain a value for Cq in the case of random periodic functions given by 
random Fourier series with independent Gaussian coefficients (Theorem 2.7) 
and indicate the implications for random trigonometric polynomials. 

The corresponding result for two-dimensional domains. Theorem 3.8, is 
considerably more subtle to state and left to Section 3. In particular, as 
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is explained in Section 3.1, the obvious generalization of the assumption 
in (2) leads to a suboptimal choice of exponent for M. We remark that the 
results of Theorem 3.8 depend on two constants Ci and C2 associated with 
the probabilistic behavior of the function at points on the boundary and the 
interior of G, respectively. Values for these constants are computed in the 
case of random doubly periodic functions and applied to the case of random 
trigonometric polynomials in two variables. 

Not surprisingly, our results require several technical probabilistic esti- 
mates which are presented in Section 4. 

2. The one-dimensional case. This section contains the proof of Theo- 
rem 1.4, the first step of which is to provide a finite criterion against which 
one can determine if the homology computation is correct. Using this cri- 
terion, a value for Cq in the case of random periodic functions given by 
random Fourier series with independent Gaussian coefficients is computed. 
This result is then applied to the simpler setting of random trigonometric 
polynomials. The latter theoretical results are then compared against data 
obtained from numerical simulations. 

2.1. Double crossovers and homology validation. The first step in ad- 
dressing the correctness of homology computations for random functions 
is to formulate a validation criterion which is amenable to a probabilistic 
study. This will be accomplished in the current subsection, in a deterministic 
setting. 

Let G = [a,b] C M denote a compact interval and let u:G ^M. denote an 
arbitrary continuous function. We are interested in determining the homol- 
ogy of the nodal domains 

iV± = {xG [a,b]:±u{x)>0}, 

that is, we would like to determine the number of components of each nodal 
domain. As mentioned in the introduction, rather than working with A^^ di- 
rectly, we choose a discretization size M and study the associated cubical ap- 
proximations Q^f introduced in Definition 1.1. These sets depend only on the 
function values of u at the M + 1 discretization points Xk = a + k{h — a)/M, 
where k = 0, . . . ,M. For sufficiently large M and under suitable regularity 
and nondegeneracy conditions on u, one can then expect the number of 
components of Q^j and A^^ to agree. 

In the one-dimensional situation, of course, it is straightforward to charac- 
terize when the homologies of Q^^ coincide with the homologies of the nodal 
domains A^^. This is the case if, on each subinterval of the form 

• the function u has no zero at all, or alternatively, 
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• it has exactly one zero in the interior of the interval and different signs 
on either side of the zero. 

In other words, if on each interval (xi^jX^+i) between two consecutive dis- 
cretization points, the function u has at most one zero and if u is nonzero 
at all Xk, then the homology computation is correct. This shows that in the 
context of random functions, we have to find sharp upper bounds on the 
probability of having at least two zeros (counting multiplicities) in a given 
interval. One would expect that for small intervals, this probability will be 
small. 

Unfortunately, the above characterization in terms of zeros is not directly 
applicable to random functions since in this context, the condition has to 
be formulated in terms of the finite-dimensional distributions of the random 
function. One possible approach is based on the following concept. 

Definition 2.1 (Double crossover). Let u:[a,b] ^ M be a continuous 
function and let [a, (5] C [a, b] be arbitrary. Then u has a double crossover 
on [a, 13] if 



for one choice of the sign o" E {±1}. 

Since the existence of a double crossover only depends on three function 
values of the function u, the probability for the occurrence of a double 
crossover can be estimated for random functions. At first glance, however, 
the obvious connection between double crossovers and zeros is of little use 
for our problem. Since the existence of a double crossover on [a,/3] implies 
the existence of at least two zeros, the double crossover probability can 
only provide a lower bound on the probability of having at least two zeros 
in [a,P] — yet we are interested in an upper bound. In order to accomplish 
this, we use a device which was introduced in [15]. 

Definition 2.2 {Dyadic subintervals, admissible intervals). Let u: [a,b] - 
M be a continuous function and let J = [a,/3] C [a,b] be arbitrary. 

• The dyadic points in the interval J = [a,/3] are defined as 

k 

dnk = Oi + iP - a) for all A; = 0, . . . , 2" and n G Nq. 

2" 

The dyadic subintervals of J are the intervals [dn,k,dn,k+i] for all k = 
0,...,2''-l and nGNo. 

• The interval J = [a, /3] is called admissible for u if the function u does not 
have a double crossover on any of the dyadic subintervals of J. 



(4) 





and a ■ u{p) > 
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The notion of admissibility depends only on the signs of the function 
values at the countable set of dyadic points in the interval, which makes 
it ideal for a probabilistic setting. Moreover, it allows us to determine the 
locations of zeros of a continuous function. We begin by presenting a criterion 
which rules out sign changes. 

Lemma 2.3. Let u:J be a continuous function on J = [a,/3] and 
assume that J is admissible for u. Furthermore, suppose that both u{a) > 
and u{P) > 0. Then for all x G J, we have u{x) > 0. Analogously, the 
inequalities u{a) < and u{[3) < imply that u{x) <0 on J . 

Proof. Assume that u{a) > and u{P) > 0. Since the function u is 
continuous, it suffices to verify inductively that u{dn,k) > at all dyadic 
points dn,k- According to our assumption, we have u{do^k) > for = 0, 1. 

Now, assume that for some n G No, we have shown u{dn,k) > for all k = 
0, . . . , 2". Since J is admissible, the function u has no double crossover on the 
dyadic subinterval [dn,k,dn,k+i] for any /c = 0, . . . , 2" — 1. With u{dn,k) > 
and u{dn,k+i) > 0, this immediately implies that u{{dn,k + '^n,/c+i)/2) > 
and this in turn implies that > for all £ = 0, ... , 2"^"^. □ 

The next lemma addresses the possibility of a sign change on an admissible 
interval. 



Lemma 2.4. Let tx: J— >M be a continuous function on J = [a,/3] and 
assume that J is admissible for u. Furthermore, suppose that u{a) > > 
u{(3). Then there exists a point x* G such that 

u{x) > for all a < X < X* and u{x) < for all x* <x< (3. 

An analogous result holds if u{a) < < u{P). 

Proof. We recursively introduce two sequences {ak)k<^No and (/3fc)fcGNo 
in J as follows. Let ao = a, (3^ = (3 and for A; G No, set 

A a Oik + /3k ., ((^k + l3k\.f. 
ak+i = ak and = ^ if u( ^ I < 0, 

ttfe+i = t; and I3k+i = f3k, ii u[ > 0. 



Then the sequence (ok) is increasing and the sequence {f3k) is decreasing. 
In addition, we have both Pk — o^k = {P — a)/2^ and u{ak) > > u{Pk) for 
all A: G Nq. Using Lemma 2.3, it is also straightforward to verify that for all 
k G No, we have both u{x) > on [a,afc] and u{x) < on [Pk,P]. The result 
now follows with x* = linik-^oo o^k- ^ 
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Note that in the situation of the above lemma, the function u can certainly 
have more than one zero. In addition to the zero x*, it is possible that u 
has zeros in both (a,a;*) and provided the function does not change 

sign at these zeros. In other words, any additional zeros have to be double 
zeros. 

By combining the above two lemmas, we can now formulate our validation 
criterion for homology computations in one dimension. 

Proposition 2.5 {Validation criterion). Let it:/— >M he a continu- 
ous function on the compact interval I = [a, h] and let M £N be arbitrary. 
Let denote the nodal domains of u and let denote their cubical ap- 
proximations, as in Definition 1.1. Furthermore, assume that the following 
hold: 

(a) the function u is nonzero at all grid points for k = 0, . . . , M ; 

(b) the function u has no double zero in (a, 6), that is, if x £ (a, 6) is a 
zero of u, then u attains both positive and negative function values in every 
neighborhood of x; 

(c) for every k = 0, . . . , M — 1, the interval [xk,Xk+i] between consecutive 
discretization points is admissible for u in the sense of Definition 2.2. 

We then have 

that is, the homologies of the nodal domains and of their cubical approxima- 
tions coincide. 

Proof. Consider the interval [x^, x^+i]. If n(xfc) and u{xf^^i) are both 
positive, then by Lemma 2.3 and (b), the function u must be strictly posi- 
tive on [xfc, Similarly, negative function values force u to be negative 
on [xk,Xk+i]. Finally, if the function values at the endpoints have opposite 
signs, then Lemma 2.4 and (b) imply the existence of a unique zero x* in 
{xk,Xk+i). From this, the result follows easily. □ 

Proposition 2.5 leads to a straightforward proof of Theorem 1.4. 

Proof of Theorem 1.4. Fix a discretization size M G N and for 
^ G {0, . . . , M — 1}, consider the interval Ii = [xi,Xi+i] between two con- 
secutive discretization points, as defined in Definition 1.1. The interval Ii is 
not admissible if the function u has a double crossover on one of its dyadic 
subintervals. If we now denote the dyadic points in li by dn,k, as in Defini- 
tion 2.2, then together with (2) and 6 = {b — a)/M , we obtain the estimate 

P{/£ is not admissible} < ^ ^ P{n has a double crossover in [dn,k,dn,k+i]} 

n=0 k=0 
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oo 2"-l 

^ E E (p+iK,fc,V2")+P-i«fe,V2")) 

n=0 k=0 
n=0 A;=0 

According to Proposition 2.5 and Assumptions (Al) and (A2), the probabil- 
ity that the homologies of N^{u;) and Q^{lv) differ can then be estimated 
as 

M 

1 - F{H^{N^) ^ H^iQfj)} < F{0 is a critical value} + ^P{ti(xf) = 0} 

e=Q 

+ ^ ^{1^ is not admissible} 

3 



^ 8Co(6-a 
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which completes the proof of the theorem. □ 

Note that the above proof could easily be adapted to also cover nonuni- 
form discretizations of the domain G. 



2.2. Random periodic functions. In order to apply Theorem 1.4 to spe- 
cific random sums, we must verify the assumption in (2). For the case of 
periodic random functions, this is demonstrated in the following. 

Assumption 2.6. Let ($7,J^,P) denote a probability space, define G = 
[0, L] for some fixed L > and consider the random Fourier series u:G xVt ^ 
M defined by 

oo 

(5) u{x,uj) = ^ak- 
fc=o 

We assume that the following hold: 

(a) the Gaussian random variables : — > M, ^ E No, in (5) are defined 
over the common probability space (O,^, P) and are independent and nor- 
mally distributed with mean and variance 1; 

(b) the constants a^. A; E No, in (5) are arbitrary real numbers such that 
at least two of them do not vanish and such that 

oo 

k^a\ < oo. 

fe=o 



92k{i^) ■ cos 



2irkx 
L 



+ g2fc-i(w) • sin 



2TTkx\ 
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We would like to point out that in the above assumption, the constants 
are directly related to smoothness properties of the random function u. More 
precisely, one can show that if 

oo 

k'^^af. < oo for some p> 0, 

k=0 

then P-almost all realizations u{-,uj) are contained in the Holder space C"'[0, L] 
for any real < q < p; see, for example, [22], Section 7.4. Moreover, if we 
define 

oo 

(6) A, = J2k''al 

k=0 

then it is straightforward to verify that, in fact, 

where E denotes the expected value of a random variable over (17, J^, P). In 
other words, the constant contains averaged information on the L^{0, L)- 
norm of the £th derivative of the random function u. 

Random periodic functions as in (5) are particularly convenient to study 
due to special properties of their spatial covariance functions. To see this, 
consider the spatial covariance function of (5), which is given by 

^ 2 27rk{x-y) 
2^«fc-cos . 



(7) R{x,y) = Eu{x)u{y) = 2_^ak- cos 

k=0 

In other words, under the above assumptions, the random function defined 
in (5) is a homogeneous random field in the sense of [1] — this simplifies some 
of the considerations significantly. 

Now consider the problem of determining the homology of the nodal do- 
mains of u using a spatial discretization of size M . 

Theorem 2.7 (Random periodic functions in ID). We consider the ran- 
dom L-periodic function u : G x — > R defined in (5), where G = [0, L] with 
L> 0, and suppose that Assumption 2.6 holds. 

Let M denote an arbitrary positive integer and, as in Definition 1.1, 
let Q^{ijj) denote the cubical approximations of the random nodal domains 
N^{lo) of u{-,uj) which are used for the computation of the homologies of 
N^{lo). Then the probability that the homology of the random nodal domains 
N^{uj) is computed correctly with the discretization of size M satisfies 



(8) : H.iQ-)} > 1 - ^ . ^||^ ^0{^). 

where the constants A^ are defined as in (6). 
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Proof. According to our assumptions, the random variable u{x, •) : — > 
M is normally distributed with mean and variance J^'k^o "^fc 7^ each 
X £ G, which immediately implies (Al). Furthermore, (A2) follows readily 
from [1], Theorem 3.2.1. Thus in order to apply Theorem 1.4, we need only 
verify (2). 

We employ Proposition 4.1 with n = 3 and sign vector (si, S2, S3) = (1, —1, 1). 
Let X £ G he arbitrary, but fixed, and consider the 5-dependent three- 
dimensional random vector T{5) = {Ti{6),T2{5),T3{6)) -.n defined by 



Ti{6)=u{x), T2i6)=u[x + 



and ^3(5) 



Then T{6) is a Gaussian random variable with mean G 
matrix C{6) given by 

" r(0) r{6/2) r{S) ' 
r{5/2) r(0) r{6/2) 
_ r{5) r{5/2) r(0) 

where, in view of (7), we use the abbreviation 



u[x + 5). 

^ and covariance 



C{5) 



r{5) = R{y + 5, y) = ^ a|. • cos 

fc=0 



27rM 
L 



For 5 — > 0, the even function r can be expanded as 

r{5) = RQ-?^.5^ + ?^.5^-^.5^ + 0{5^) with i?, : 
^ ^ 2 24 720 ^ ' 

and Assumption 2.6(b) implies that 



2f 



27ry 



(9) 



i?o/0, i?i/0 and RoR2-Ri>0. 



For the last inequality, one need only observe that the Cauchy-Schwarz 
inequality implies i2f < R0R2, and that equality can easily be excluded since 
two of the Ofc are nonzero. Using the above expansion, the determinant of 
the covariance matrix can be written as 

(10) det C{6) = ^ ■ {R0R2 - Ri) ■ 5^ + 0{5^), 

64 

that is, the covariance matrix is positive definite for sufficiently small 5 > 
and Proposition 4.1(i) is satisfied. Furthermore, by applying the Newton 
polygon method [27, 29] to the characteristic polynomial det(C((5) — AI), 
it can be shown that in the limit (5 — > 0, the three eigenvalues Afc((5), for 
A; = l,2,3, of C{5) satisfy 



Ai(<5) 



i?oi?2 



^1 



96i?o 
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(11) 



A2W 



Ri 
2 
3i?o 



■6^ + 0(6^) 
0{5). 



This establishes the vahdity of Proposition 4.1(iii). In order to determine 
the asymptotic behavior of the eigenvector corresponding to Ai, we consider 
the adjoint of the covariance matrix, whose expansion is given by 



adjC(5) 



RqRi 



1 



1 



1 



6"^ + 0(6^). 



The constant coefficient matrix has the double eigenvalue 0, as well as the 
positive eigenvalue 6 with associated unnormalized eigenvector (1,-2,1)*. 
Since the eigenspace for the largest eigenvalue of the adjoint matrix coincides 
with the eigenspace for the eigenvalue Xi{d) of C{5), the simplicity of these 
eigenvalues shows that we can choose an eigenvector vi{5) for Xi{5) which 
satisfies vi{d) — > (1, — 2, l)*/6^/^ for (5 — > and therefore Proposition 4.1(ii) 
is satisfied. Applying Proposition 4.1, we finally obtain 



lim p+i{x,6) 



I det C{5) _ sVg 



Ai(5)3 Svr ■ 

In combination with (9), (10) and (11), this limit implies that 



Pa{x,6) 



UStt 



Rq Ri 



6^ + 0{6^ 



16L3 



43/2 .1/2 



■ + 0{6' 



for a = 1, the analogous statement for a = —1 following from a symmetry ar- 
gument. Thus (2) is also satisfied and Theorem 2.7 follows from our abstract 
result, Theorem 1.4. □ 



Very much in the spirit of [25], our result furnishes an explicit bound on 
the likelihood of computing the correct homology. The bound depends on 
the discretization size M and global properties of the underlying function u. 
Unlike [25], however, the necessary information on u can be easily computed. 
In fact, it is given explicitly in terms of the coefficient sequence (a^) which, 
in turn, is related to smoothness properties of u. Note also that due to the 
asymptotic nature of our results, the estimate in (8) contains an additional 
0(1/M3) term. While in principle it is also possible to derive explicit bounds 
for this term, the formulas involved are rather complicated. In practice, 
however, the term can already be neglected for relatively small discretization 
sizes. For more details, we refer the reader to [11, 12]. 
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2.3. Random trigonometric polynomials. We now demonstrate the impli- 
cations of our main results for the case of random trigonometric polynomials. 
Let > 3 denote an arbitrary integer and let L = 2tt. If, in the situation of 
Theorem 2.7, we choose the coefficients specifically as = 1 for 1 < /c < A^, 
ao = and = for k > N, then (5) represents a random trigonometric 
polynomial, 

TV 

(12) u{x,Lo) = ^(st2fc(w) • cos(A;j;) + g2k-ii^^) ■ sm{kx)). 

k=l 

The asymptotic behavior of the number of zeros of such functions has been 
studied extensively and therefore provides us with a means for testing The- 
orem 2.7. Specifically, one can show that as — > oo, most trigonometric 
polynomials of the form (12) have on the order of 2N/3^^'^ real zeros in the 
interval [0,27r]; see, for example, [4, 16, 17]. Thus in order for our homology 
computation to be correct, we must at least ensure that the discretization 
size M satisfies 1/M = 0{1/N). Due to the almost certain occurrence of 
zeros which are more closely spaced, it seems implausible that this bound 
would be optimal. Now, let be defined as in (6). Then one can easily 
show that 



.3/2 ,1/2 - - '^^^^ + + l)(2Ar + 1) ~ ^ . AT^ 

This implies that in order for the homology computation to be accurate with 
high confidence, we must choose the discretization size M in such a way that 

M~Af^/2 forAf^oo. 

We conjecture that this proportionality is the correct asymptotic result, as is 
supported by the numerical results shown in Figure 4. For various values of A^ 
between 5 and 1000, we computed three different quantities: the expected 
number of zeros (lowermost curve), the expected minimal distance between 
two consecutive zeros of the random polynomial (second curve from bottom) 
and the integer M for which 95% of the generated random polynomials have 
a minimal distance between consecutive zeros which is larger that 2tt/M 
(uppermost curve). For comparison, the prediction of Theorem 2.7 which 
would ensure a correctness probability of at least 95% is given by the curve, 
which lies exactly below the uppermost curve and shows the same asymptotic 
growth behavior. 

Numerically validated computational results which also treat the case of 
linear and nonlinear evolution equations (such as the Cahn-Hilliard models 
described in the introduction) will be presented in [12]. 
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3. The two-dimensional case. In this section, we address the case of two- 
dimensional base domains G cM? and are interested in the topology of 
the nodal domains of a function u:G^M. For simplicity, it is assumed 
that G is a square domain, even though more general domains could also be 
treated. As in the one-dimensional situation, we present both general results 
and specific applications, in particular, the case of doubly periodic random 
functions. 

3.1. Double crossovers and B -admissibility. In contrast to the one-dimen- 
sional situation where a nodal domain of a function consists of a disjoint 
union of intervals, the geometry of nodal domains can be rather complicated 
in two space dimensions. It is therefore not surprising that the method of 
Section 2 cannot be generalized directly to the higher-dimensional case. In 
this first subsection, we show that a fairly straightforward adaptation of 
the one-dimensional ideas can be used to discuss the accuracy of homology 
computations. Unfortunately, however, this approach alone will provide sub- 
optimal results and must therefore be extended significantly in the following 
subsection. 

In the one-dimensional situation, the basic ingredients of our approach 
can be summarized as follows. After recognizing the fundamental obstruc- 
tion to a valid homology computation as the occurrence of two zeros in an 
interval between gridpoints, this obstruction was reformulated in terms of 
the concept of a double crossover. The latter is a specific function value sign 
pattern that can be observed at the endpoints and midpoint of a given inter- 
val J. Moreover, by excluding this sign pattern from all dyadic subintervals 
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Fig. 5. Examples of problem situations in two dimensions. In each of the three columns, 
the function values at the grid points shown have identical sign configurations, but the 
underlying topology of the nodal domains differs. The nodal lines are shown as solid black 
lines. 

of J, we were able to introduce the notion of admissibility of J. Admis- 
sibility guarantees that the topological information for the nodal domains 
in J can be determined from the function values of the given function at the 
endpoints of J. If the function values at the endpoints have the same sign, 
then the function has no zero in J; if the signs differ, then the function has 
exactly one zero in J. 

Based on the above discussion, our goal is the identification of certain sign 
patterns of function values that could indicate mistakes in the homology 
computation. As the fundamental building block, we consider a square with 
side length 6 and, analogously to the one-dimensional case, we consider the 
function values of the underlying function u at the corners as well as at the 
midpoints of the sides and at the center of the square. This will easily enable 
us to consider dyadic subsquares as before. 

Some representative examples of problematic sign patterns are shown in 
Figure 5. Within each column, the two examples exhibit the same sign pat- 
tern, but the underlying topology of the nodal domains differs. The exam- 
ples in the first two columns are, in some sense, reminiscent of the one- 
dimensional situation since in each of these cases, we can find a double 
crossover on at least one vertical or horizontal line of grid points. This is 
not true for the examples in the last column of Figure 5. Here, there are no 
double crossovers on vertical or horizontal lines, yet the underlying topology 
still differs. On the level of the function values of u, this is due to the fact 
that, at the corners, the signs alternate. In other words, there are two double 
crossovers around the perimeter of the square. Altogether, we see that each 
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Fig. 6. For the notion of B -admissibility, the seven sign configurations shown above 
are forbidden. Each image corresponds to two sign patterns: one m which the dark points 
correspond to positive function values and the light points correspond to negative ones, and 
one for the opposite situation. 



of the sign configurations shown in Figure 5 contains one of the seven sign 
patterns depicted in Figure 6. 

It turns out that excluding the patterns in Figure 6 is sufficient to extend 
the general procedure of Section 2. For this, we need the following definitions. 



Definition 3.1 {Dyadic subsquares). Let J := [a, a + S]x[(3, (3 + 5] C 
and denote the dyadic points in the interval [0,6] in the sense of Defini- 
tion 2.2, by dn,k- Then the dyadic points in the square J are the points 

dn,k,i = {a + dn,k, l3 + dn/) G for all fc, £ = 0, . . . , 2" and n G Nq. 

The dyadic subsquares of J are the squares dn,k,e + [0,(^/2"]^ for all k,£ = 
0,...,2" - 1 and nG No. 

Definition 3.2 {B- admissibility). Let urG— >M denote an arbitrary 
function and let J C G denote a square, as in Definition 3.1. Then J is 
called B- admissible for u, if none of its dyadic subsquares contain any of 
the function value sign configurations shown in Figure 6. 

Using the concept of B-admissibility, one can validate homology more or 
less analogously to the one-dimensional case. To demonstrate this, consider 
Figure 7. If we assume that a given square J of side length 6 > (shown in 
pink) is i?-admissible for u, then the figure shows all possible function value 
sign configurations which can be observed at the dyadic points di^k,£> where 
k,i = 0, 1, 2, up to symmetry actions of the square. These configurations are 
grouped according to the sign configurations at the corners of the square. 
We now make two important observations: 
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Fig. 7. Possible sign patterns for the function values of u at the dyadic points di,fe,f, 
where k, I = 0,1,2, if the underlying square is B -admissible. In addition to the patterns 
shown, all sign patterns which can be obtained by applying an element of the symmetry 
group of the square can also occur, resulting in a total of 66 possible patterns. 

• if all function values of u at the four corners of J have the same sign, then 
all function values at the nine depicted dyadic points must necessarily 
have the same sign; 

• if both signs can be observed at the corners of J, then there are exactly 
two sides of the square with contain both positive and negative function 
values. 

The first observation is reminiscent of Lemma 2.3. It shows that if J is 
i?-admissible for u and u is positive at the corners of J, then u cannot take 
negative function values on J. Thus the first observation can be used to 
identify subsets of one or the other nodal domain. In combination with the 
second observation, it therefore seems plausible to expect that the nodal 
line in a i?-admissible square can be pinned down. In fact, under suitable 
nondegeneracy assumptions on n, it is possible to prove the following result. 

Proposition 3.3. Let G C denote a square, letu:G beC'^ such 
that is not a critical value of u and let J C G be a B-admissible square 
for u. Then the following hold. 

(a) If u is strictly positive at the corners of J , then u is strictly positive 
on J. Similarly, if u is strictly negative at the corners of J , then u is strictly 
negative on J. 

(b) If u takes both positive and negative function values at the corners 
of J, then the nodal line of u in J is a simple smooth curve which connects 
one side of J with another side of J . 
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Now, assume the situation of Definition 1.1 with a positive integer M. If 
all subsquares which are generated by adjacent gridpoints of the equidis- 
tant M -discretization of G are B-admissible for u, then we have 

if,(iV±)^if,(Q±), 

that is, the homologies of the nodal domains and of their cubical approx- 
imations Q^j coincide. 

The proof of the above result can easily be formulated using ideas similar 
to those for the one-dimensional case and is therefore omitted. 

At first glance, Proposition 3.3 is exactly the validation criterion which 
allows us to consider the two-dimensional random field case. It is, in fact, 
possible to proceed as in Section 2 and formulate a probabilistic estimate 
based on the above proposition. Unfortunately, the resulting theorem pro- 
vides highly suboptimal probability estimates. This suboptimality stems 
from the fact that the three-point sign configurations in Figure 6 are not 
sufficiently unlikely. It turns out that in many applications, each of these con- 
figurations has probability of order 0{6^), where 6 denotes the side length 
of the square J. On the other hand, the four-point configuration in Figure 6 
has probability of order 0((5^). Together, the probability for ^-admissibility 
is of order 0{5^). For the homology validation, we must ensure that each of 
the subsquares of G are i?-admissible and that (5 ~ 1/M. This finally 
gives a probability estimate of order 0(1/M), in contrast to the 0(1/M^) 
estimate in one dimension. 

3.2. Excluded sign patterns and I -admissibility. The discussion of the 
last paragraph indicates two things. On one hand, one has to avoid the 
one-dimensional double crossover concept in higher dimensions since its oc- 
currence is too likely. On the other hand, the four-point configuration of 
Figure 6 does exhibit the correct probability estimate. If we were able to 
define a notion of admissibility which has probability of order 0{6^), then 
the resulting homology validation result would exhibit the desired quadratic 
order in M. 

As a first step toward a suitable admissibility condition, we return to 
the four-point configuration of Figure 6. It was mentioned above that in 
many applications, the probability of this sign pattern turns out to be of 
order 0(6'^) — one would expect that nondegenerate affine transformations of 
this pattern are similarly rare. This leads naturally to the following notion. 

Definition 3.4 {^-admissibility). Let u:G^R and J = [a,a + 5] x 
[P, P-\-5] C G. Then J is called I 4^- admissible for u if it does not contain any 
of the function value sign configurations shown in Figure 8. 
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Fig. 8. For the notion of li-admissibility, the sixteen sign configurations shown above 
are forbidden. As in Figure 6, each image corresponds to two sign patterns. 



We are still interested in a recursive approach to homology validation, as 
described in Sections 2 and 3.1. It is therefore natural to determine which 
of the 512 sign patterns on the dyadic points di^k/, where k,i = 0, 1,2, can 
occur if the underlying square is /4-admissible. An exhaustive search shows 
that only 92 patterns remain possible. Up to symmetry actions of the square, 
these configurations are shown in Figure 9. 

In contrast to the situation described in Section 3.1, not all patterns that 
are possible under /4-admissibility permit homology validation. The prime 
example can be seen in the first row of Figure 9. Even though the signs at 
the corners of the square are all the same, after one level of refinement, the 
opposite sign appears at the center of the square. This example shows that 
the concept of /4-admissibility alone cannot be used for homology validation, 
that is, we must exclude more subconfigurations. At the same time, we 
must ensure that the probability estimates for the additionally excluded 
configurations are sufficiently small. As we will see later, it can be shown 
that in the context of random periodic functions, the last sign configuration 
in the first row of Figure 9 exhibits an 0{6^)-hound. In fact, we need not 
consider the full nine-point pattern, but can rather consider the five-point 
subpattern shown in Figure 10. This leads to the following definition. 
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Fig. 9. Possible sign patterns for the function values of u at the dyadic points di^k,i, 
where k, I = 0,1,2, if the underlying square is 14-admissible. In addition to the patterns 
shown, all sign patterns which can be obtained by applying an element of the symmetry 
group of the square can also occur, resulting in a total of 92 possible patterns. 

Definition 3.5 admissibility). Let u:G^R and J = [a,a + 5] x 
[(3, (3 + 5] C G. Then J is called I^-admissible for u if it does not contain the 
function value sign configuration shown in Figure 10. 

While combining the notions of I4- and /5-admissibility does bring us 
closer to homology validation, problematic configurations can still be identi- 
fied. This is demonstrated in Figure 11 by considering groups of two adjacent 
/4- and /5-admissible squares. In each case, the sign patterns at the corners 
of the original squares do not suffice to determine the topology after one 
refinement level. 

The topological changes shown in Figure 11 are introduced by specific 
patterns in Figure 9, namely the second nine-point pattern in the first row 
and the first nine-point pattern in the last row. Again in the context of 
random periodic functions, one can determine their likelihood, which pro- 
vides only the insufficient 0(5'^)-bound. It is therefore impossible to exclude 
these two patterns directly, or any of their subpatterns. Instead, we follow 
a different strategy. Note that the two rightmost patterns in Figure 11 can 
be avoided if the square with center at the midpoint of the common edge 




Fig. 10. For the notion of -admissibility, the five-point sign configuration shown above 
is forbidden. As in Figure 6, this image corresponds to two sign patterns. 
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Fig. 11. Sample configurations which cannot be detected by using both Ia- 
Is-admissibility, but which exhibit different topological properties. 



and 



between the two original squares is both /4- and /5-admissible. Combining 
this with the concept of dyadic subsquares (as in Section 3.1) leads to the 
following, final, definition. 

Definition 3.6 {I -admissibility). Consider u:G^M and J=[a,a + 
5] X [/?,/? + (5] such that the closed 5/2-neighborhood of J is contained in G. 
Then the square J is called I-admissible for u if for every dyadic sub- 
square J* of J with side length 6* = (5/2", n G No, the following five squares 
are both /4-admissible and /5-admissible for u: 

• the dyadic subsquare J* itself; 

• the four shifted squares which are obtained by translating J* horizontally 
or vertically by 6* /2 in either direction, as shown in Figure 12. 

As mentioned before, it will be shown later that in many situations, it 
is possible to derive an 0((5^)-bound for /-admissibility. Note, however, 
that this improved bound comes at an additional cost. The notion of /- 
admissibility is not determined solely from the function values of u at the 
dyadic points in J. Rather, one must also include knowledge from certain 
dyadic points in neighboring squares and this renders /-admissibility im- 
practical at the boundary of G. In order to resolve this final issue, we must 
combine both B- and /-admissibility in a suitable way. This is presented in 
the next subsection. 



3.3. A criterion for homology validation. In the current subsection, the 
notions of B- and /-admissibility are combined to formulate and prove a 
homology validation criterion in two space dimensions which is amenable to 
a probabilistic treatment. 

While Theorem 1.2 provides a criterion under which it can be guaranteed 
that for sufficiently large M, the homology of the nodal domains agrees with 
that of the cubical approximation, it does not provide insight into how we 
would determine a minimum value of M. For this, we must consider the 
notions of B- and /-admissibility. 
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Proposition 3.7 (Validation criterion). Let G C denote a square do- 
main with sides parallel to the coordinate axes and let u:G ^M. be a twice 
continuously differentiate function for which is a regular value. Let 
denote the nodal domains of u and let Q^j denote their cubical approxima- 
tions (as in Definition 1.1) for some fixed M £ N. In addition, suppose that 
for each of the closed subsquares J formed by adjacent gridpoints in G, 
u is nonzero at all dyadic subpoints, and that one of the following holds: 

(a) if the square J lies at the boundary ofG, then J is B -admissible for u; 

(b) if the square J lies in the interior of G, then J is I-admissible for u. 

We then have 



that is, the homologies of the nodal domains and of their cubical approxima- 
tions coincide. 

Proof. Let M be fixed as in the formulation of the proposition. In view 
of Theorem 1.2, it suffices to show that for all values of n G N, the homologies 
of Q^.i.2n and coincide since then we need only choose n sufficiently large 
for the lemma to hold. In order to establish this invariance of the homologies 
inductively, we only consider the homologies of Q2M Q%i- Subsequent 
inductive steps can be carried out completely analogously. 

Now consider, without loss of generality, a connected component K of the 
set We can associate with K a collection of positive grid points P in 



Fig. 12. For the notion of I -admissibility of a square J, we require that every dyadic 
subsquare J* (shown with shaded background) of J, as well as the four shifted versions 
of J* shown on the right, are both I4,- and I^-admissible, In the figure, the location of the 
shifted squares is indicated by the five-point configuration of Figure 10. Nevertheless, all 
configurations in Figure 8 must be excluded from the five indicated squares. 



H,{N^)^H,{Q%), 






26 



K. MISCHAIKOW AND T. WANNER 




Fig. 13. Possible sign patterns for the function values of u at the dyadic points di^k,i, 
where k,£ = 0,1,2, if the underlying sguare is I -admissible. In addition to the patterns 
shown, all sign patterns which can be obtained by applying an element of the symmetry 
group of the square can also occur, resulting in a total of 90 possible patterns. However, the 
two configurations with patterned backgrounds can only occur in certain situations which 
are described in more detail in Figure 14- 

the underlying discretization of G, that is, grid points at which the func- 
tion u takes positive function values. We must show that the insertion of 
intermediate grid points in the underlying grid does not change the topology 
of the component K in the following sense: The corresponding component 
in Q2M merged with another component and the number of holes 

in K has remained constant. 

In order to perform the step from the given discretization of size M to the 
one of size 2M, we must establish which nine-point sign configurations are 
possible for every given four-point sign configuration at the corners of a dis- 
cretization square. For the concept of i3-admissibility, this has already been 
accomplished in Figure 7. The resulting patterns immediately show that for 
a i?-admissible square, passing from M to 2M cannot change the topological 
properties since no new connections are formed and positive/negative grid 
points are only introduced adjacent to positive/negative components. 

As for the concept of /-admissibility, we must be more careful. In Fig- 
ure 13, we collect all possible nine-point sign configurations which could be 
observed at the next refinement level. With the exception of the two con- 
figurations with patterned backgrounds, the same comments as in the last 
paragraph apply to these patterns and they therefore cannot introduce a 
change in topology. 

Now, consider the two remaining nine-point configurations in Figure 13, 
that is, the configurations with patterned backgrounds. Recall that our defi- 
nition of /-admissibility includes information about neighboring squares and 
this implies that these two configurations can only occur if at the edge of 
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Fig. 14. The two configurations in Figure 13 with patterned backgrounds can only occur 
if the neighboring square exhibits a specific sign configuration. For the pattern in the first 
row of Figure 13, this is shown in the left-hand diagram; for the configuration in the third 
row see the right-hand diagram. 

the square with the newly generated double crossover, the adjacent square 
exhibits the specific sign pattern shown in Figure 14. Thus the newly cre- 
ated double crossovers can only occur at the boundary of the component K 
and if they appear, they do not create any new connections or introduce 
any change in topology. This fact also implies that in the interior of no 
negative gridpoints can be introduced. In addition, note that if the adjacent 
square is, in fact, a S-admissible square, then the two patterned configura- 
tions cannot occur at all due to the generation of a double crossover. These 
facts imply that the topology of the component K will not be changed by 
passing from Q2M proof of the proposition is complete. □ 

We are finally in a position to present an abstract version of a probabilistic 
homology validation result for two-dimensional random fields u on a square 
domain G CM? . 

Theorem 3.8 (Abstract 2D estimate). Consider a probability space (17, J^, 
P) and the square domain G = [0, L]^ C M^. Let « : G x O — > M denote a ran- 
dom field over (fi,^, P) satisfying Assumptions (Al) and (A2) such that for 
¥-almost all lo the function u{-,uj) : G — > M is twice continuously differ- 
entiable. For each w G f], denote the nodal domains of u{-,uj) by N'^{u)cG 
and denote their cubical approximations (as in Definition 1.1) by Q^j{uj). 

For X = {xi,X2) and 6 > such that J = [xi,xi -\- 6] x [x2,X2 + 6] C G, 
consider the following events in T : 

(i) let Eb{x,6) denote the set of all lo £ Q for which u{-,uj) exhibits at 
least one of the seven sign patterns in Figure 6 at the nine dyadic points dn^^ 
/or £,m = 0,1,2; 

(ii) let Ei{x,6) denote the set of all uj £ for which u{-,uj) exhibits at 
least one of the seventeen sign patterns in Figures 8 and 10 at the nine 
dyadic points di^i^m for i,m = 0,1,2. 
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Assume that there exist positive constants Ci and C2 such that 

(13) F{Eb{x,6))<Ci-6''^ and F{Ei{x,6)) < C2 ■ 6^ 

for all X £ G and 5 > for which J lies in G. 

Then for every discretization size M , the probability that the homologies 
of N^iuj) and Q^j{lo) coincide satisfies 



(14) F{H,{N^)^H,iQ^,)}>l 



Mil > i 3^^2 



Proof. Fix a discretization size M S N. Then the discretization points 
subdivide the square G into subsquares, each with side length 5 = L/M. 
Let J denote one of these subsquares which hes at the boundary. In order 
to estimate the probabihty of J not being i?-admissible, we must consider 
dyadic subdivisions, as in Definition 3.1. At the nth subdivision level, there 
are 4" dyadic subsquares with sides of length 5/2". According to (13), the 
probability that a particular one of these dyadic subsquares contains one 
of the forbidden sign configurations of Figure 6 is bounded by Ci • ((5/2")^. 
Summing over all dyadic subsquares and subdivision levels n G No, we then 
obtain 



[J is not i?-admissible} < ^ ^ Ci • f—^ 

n=0 k=0 



M3 



Now, assume that J is contained in the interior of G. As above, (13) implies 
that 

(15) P{ J is not /-admissible} < ^ ^ 5C2 ' ^ ^ 

n=0 k=0 



2" J 3M4 



Note the additional factor of 5 due to the fact that for each of the dyadic 
subsquares, we must also allow for the excluded sign configurations in the 
four shifted squares; see Definition 3.6. 

According to Proposition 3.7, if w G $7 is such that the homologies of N^{uj) 
and Q^,^{uj) differ, then one of the following statements must be true: 

• the function u{-,uj) vanishes at one of the gridpoints or at one of the 
dyadic subpoints of the grid subsquares; 

• the function u{-,u!) has a double zero; 

• one of the subsquares J which lies at the boundary of G is not in- 
admissible; 

• one of the interior subsquares J is not /-admissible. 
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Due to Assumptions (Al) and (A2), the first two events have probabihty 
zero and we therefore finally obtain 

1 - P{H^iN^) - H.iQ^)} < (M^ - (M - 2)2) • ^ + (M - 2f 



M3 ^ ^ 3M4 
24CiL3 + 20C2L^ 
- UP ' 

which completes the proof of the theorem. □ 

The above result could be simplified for the situation of doubly-periodic 
functions. If the domain of such a function is identified with a torus and 
we are interested in the homology of the nodal domains as subsets of this 
torus, then we can reformulate the theorem by employing /-admissibility 
alone. Furthermore, due to the periodicity of the domain, we can replace 
the additional factor of 5 in (15) by 3 since it now suffices to check each 
square, as well as its down-shift and its right-shift. Thus for doubly periodic 
functions, the final probability bound reads 

AC T ^ 

P{//.(iV±)-i/.(Q±)}>l-^, 

where C2 is as in (13). Finally, as will be demonstrated in the next subsection, 
the bounds in (13) are optimal for doubly periodic random functions. 

3.4. Random periodic functions. As in the one-dimensional case, we close 
this section by providing explicit results for the important special case of 
doubly periodic functions, that is, for classical random Fourier series. 

Assumption 3.9. We consider random Fourier series on G= [0,L]2 of 
the form 

/ 27rA:xi 2TTix2 

u[x,uj)= 2^ Ofc/- { gk,e,i[^) cos — ; — cos 



kl=0 



(16) 



+ 9k,£,2[^) cos 



L L 

2'Kkxi . 2t:Ix2 



L 


sm 


L 


2TTkxi 




21X1X2 


cos 




L 




L 


2'irkxi 
I — - — 


■ sin 


2TTtX2 



L L 



which satisfy the following: 



(a) The random variables gk/,m in (16) are defined over a probability 
space {Q,J-',¥) and are independent and normally distributed with mean 
and variance 1. 
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(b) there exist positive ki,ii G N and nonnegative k2,i2 G No which sat- 
isfy ki / k2 and 7^ £2, as weU as fcf + £f / + ^2 ; such that both akj/-^ 
and ak2/2 nonzero. In addition, suppose that 

00 

Y^ik' + i')al,<oo. 
k,e=o 



As before, the summability condition in Assumption 3.9(b) is related to 
smoothness properties of the random field u. In fact, it guarantees that 
the function u{-,uj) has continuous partial derivatives up to order two for 
P-almost all u; G 17. Moreover, if we define 

00 

(17) A,,, = Y: k^'^^'ale, 

k,e=o 

then we can readily show that 

where E denotes the expected value of a random variable over (0,JF, P). 

As in the one-dimensional case, doubly periodic functions simplify our 
discussion, in that their spatial covariance function is shift-invariant. More 
precisely, the spatial covariance function of (16) which is given by 

(18) R{x, y) = 'Eu{x)u{y) = ^ a\f^ ■ cos — ^ . cos ^^^^^^ — 

k,i=o ^ ^ 

Adopting this setting, we obtain the following result. 



Theorem 3.10 (Random doubly periodic functions in 2D). Consider 
the random Fourier series u defined in (16) and suppose that Assumption 3.9 
is satisfied. Let M denote an arbitrary positive integer and let Q^j{uj) denote 
the cubical approximations of the random nodal domains N^{lo) of u{-,lo) 
(as in Definition 1.1) which are used for the computation of the homologies 
of N^{ll)). Then the probability that the homology of the random nodal do- 
mains N^{ll!) is computed correctly with the discretization of size M satisfies 



P{//*(iV±)^F,(Q±)}>l 

(19) 



± ^ ^ 10677r2 (^2.0 + ^1,1 + ^0,2)^ 
18M2 ■ .1/2 .1/2 -1/2 .1/2 
^0,0 ^0,1 ^1,0 ^1,1 



+ 



where Ap^g was defined in (17). 
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Table 1 

Covariance matrix spectral information for the first shown sign pattern. For the second 
pattern, we must replace Roj by Re,o 



D 



det C{S) = .s^ + o{S') 

64 



Ai(<5) = 



Ro.oRo.2 — Ro.l ^4 

96J?o,o 



X3{S) = SRo,o + 0(S) 



Table 2 

Covariance matrix spectral information for the sign pattern shown 




det C((5) = Ro.oRo.iRi.oRia ■ + 0{5^) 



MS) ^^-S^ + OiS") 
, , 2_Ro.i-Ri,o 



Ro,i+Ro,i + 
Ro,i + Ro,i + 



A3 (5) = 
X4S) = 4Ro,o + 0{5)^ 



Ro,i — Ro,i 
Ro,i — Ro.i 



■5^ + 0{5'') 



vi(5) = i-(l,-l,l,-ir + 0(5) 



Proof. According to our assumptions, the random variable u{x, •) : — > 
M is normally distributed with mean and variance J2'k'e=o ^\t^^ fo'^ each 
X £ G and this implies (Al). Furthermore, (A2) follows from [1], Theo- 
rem 3.2.1. Thus in order to apply Theorem 3.8, we need only verify (13). 

As in the proof of Theorem 2.7, we employ Proposition 4.1 with n = 3 or 
n = 4 to estimate the probability of the various sign patterns, if one considers 
a square with side length 6, as in (13). Since this can basically be done as in 
the one-dimensional case, we have only collected the necessary spectral in- 
formation for the various covariance matrices in Tables 1-6. In each of these 
tables, the expansions of the determinant of the covariance matrix C{6), 
the expansions of the eigenvalues and the expansion for the eigenvector cor- 
responding to Ai((5) are given. The constants Rp^q which appear in these 
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Table 3 

Covariance matrix spectral information for the sign pattern shown 



det C{S) = ^o>°^o.i^i>°^i.i . ^8 + 0(^9) 
256 



, /c-N Ro,iRi,o 



2{Ro,i + Ro,i + \Ro,i ~ Ro,i\ 



Ro,i + Ro,i + \Ro.i — Ro,i\ ^2_|_Q^^3-j 
8 



As (5) = 

X4{S)=4Ro,o + 0(Sy 



vi(5) = i-(l,-l,l,-l)' + 0(5) 



Table 4 

Covariance matrix spectral information for the first sign pattern shown. For the secona 
pattern, we must replace Ro,i by Rifi and Rk,o by Ro,k 




detC(5) = 



Ro,iRi,oiRo,oRi,i + Ro,oR2,o — Ri,o) 
256 



A2(<5) 



MS) 



647?o,o 
Ro.iRi.o 



2Ro,i + 47ii,o + 2 y^Rl^+4Rl o 

Ro.l + 2i?i_o + Rq,! + 4_Rj Q 



Xi{S) = ARo.o + 0{S) 



t;i(5) = i-(l,-l,l,-l)*+0(5) 



formulas arise in the expansion of the spatial covariance function (18). More 
precisely, if we set 

, ^, , , 2 2TTkdi 27r^d2 

r{d) = R{y + d,y)= ^ a^ ^ ■ cos ^ • cos ■ 



kl=0 



L 



then for (i — > 0, the even function r can be expanded as 
r{d) = Ro,o-l-{Ri,odl + Ro,idl) 
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where 



R. 



■kl 



27r\2fc+2£ 

l) 



A 



ki- 



ln addition, we can easily verify that all of the leading coefficients in the 
expansions of the tables are nonzero due to Assumption 3.9(b). 

By combining the information in Tables 1 and 2 with Proposition 4.1 and 
the definition of Eb{x,6) in (13), we obtain the estimate 

r{EB{x,d)) <6 - -ttt:: ,^ , ,o + 6 ■ 



1287r 



p3/2 „l/2 
^0,0 ^0,1 



1287r 



p3/2 „l/2 
-"■0,0 ^1,0 



Table 5 

C'ovariance matrix spectral information for the first sign pattern shown. For the secona 
pattern, we must replace Ro,e by Ri,o and Rk.o by -Ro.fc 




A2(<5) 



64J?o,() 
Ro,iRi,o 



2Ro,i + Ri,o + VlMi±^L£ 

XsiS) = ■ ■ — ■ 5^ + 0{5^) 

A4(<5) =4i?o,o + 0(<5) 



«i(5) = i-(l,-l,l,-l)' + 0(5) 



Table 6 

Covariance matrix spectral information for the sign pattern shown 




R* = Ro,o{Ro,2 + 2J?i,i + i?2,o) - (i?o.i + -Ri,o) 



J?' 



SORo.o 
Ri.i 

4 



Ai(5) 

A2(5) 
A3(5) 

A4(J) 

A5(i5) = 5J?o,o + 0(<5)' 



2-Ro, 1^1.0 



-Ro,i +iii.o + 
-Ro.i +^1.0 + 



^0,1 — Ri,o 
Ro.i — Ri.o 



vi{5) {1,1,1,1 -^Y+OjS) 
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+ 2- 



127r2 ol/2„l/2„l/2 
^0,0 ^0,1 ^1,0 

-64^ ol/2 pl/2 pl/2 fr^ ' pl/2pl/2pl/2 

-"-0,0 -"-0,1 -"-1,0 -"-0,0 -"-0,1 -"-1,0 



8L3 4 1/2. 1/2. 1/2 3^4 .1/2 ,1/2 .1/2 

^0,0 ^0,1 ^1,0 ^0,0 ^0,1 ^1,0 

<4l7r253 (^2,0 + ^1,1 + ^0,2)=^/^ 

- 12L3 ■ .1/2 .1/2 ,1/2 

^0,0 ^0,1 ^1,0 

Note that the three terms after the first inequahty correspond to the bounds 
for the two patterns in Table 1 and the pattern in Table 2 (in that order); the 
additional factors 6, 6 and 2 correspond to the multiplicities of each pattern. 
For example, the first vertical pattern must be tested twice (depending on 
whether the interior point is positive or negative) on each of the three vertical 
lines. 

From the above estimate, we can readily see that the first part of (13) is 
satisfied with 

4l7r2 (^2,0 + ^1,1+^0,2)^ 



Ci 



12L3 .1/2 ,1/2 ,1/2 ,1/2 
^0,0 ^0,1 ^1,0 ^1,1 



Similarly, by combining the information in Tables 3-6 with Propositions 4.1 
and 4.2 (for the five-point sign pattern) and recalling the definition of Ej{x, 5) 
in (13), we obtain the estimate 

F(F(rS))<H ^ "^''^ I P iR2,o+Ri,ir/' 

n^nx, 0))_6- ^^2^2 • 1/2 1/2 1/2 + • 192^2 " 1/2 1/2 1/2 

-"-0,0 ^0,1 ^1,0 ^0,0 ^0,1 ^1,0 

{Ro,2 + Ri,i?/^ , , 6^ (fi2,o + 4fli.i)3/2 



+ 4 



+ 2 



1927r2 ol/2„l/2^1/2 384^2 ^l/2„l/2„l/2 

^0,0 ^0,1 -"^1,0 -"-0,0 -"-0,1 -"-1,0 

{Ro,2 + m,if/^ 

3847r2 ■ ol/2„l/2„l/2 
-"-0,0 -"-0,1 -"-1,0 

6^ {Ro,2 + 2iii,i + R2,o? 



< 



10247r2 „l/2„l/2„l/2„l/2 
-"■0,0 -"-0,1 -"-1,0 -"-1,1 

115J^ iRo,2 + Ri,i+R2,o? 

3847r2 ■ ol/2ol/2„l/2„l/2 ' 
-"-0,0 -"-0,1 -"-1,0 -"-1,1 



PROBABILISTIC VALIDATION OF HOMOLOGY COMPUTATIONS 



35 



that is, we can set 

_ IIStt^ (^2,0 + ^1,1 + ^0,2)^ 
^ 24L4 ■ .1/2 ,1/2 .1/2 .1/2 • 
^0,0 ^0,1 ^1,0 ^1,1 

Thus, (13) is also satisfied and Theorem 3.10 follows from our abstract result. 
Theorem 3.8. □ 

It is clear from the above proof that the bound in (19) is not optimal. Its 
form is chosen to point out the essential ingredients, but not to establish the 
tightest fit. We do believe, however, that it reflects the true situation fairly 
accurately. 

3.5. Random bivariate trigonometric polynomials. As in the one-dimensional 
case, we illustrate our main theorem for random doubly periodic functions 
using random trigonometric polynomials. For this, let > 3 denote an arbi- 
trary integer. If, in the situation of Theorem 3.10, we choose the coefficients 
as Uk/ = 1 ior 1 < k,i < N and a^,^ = otherwise, then (16) represents a bi- 
variate random trigonometric polynomial. Now, let be defined as in (17). 
We can then easily show that 

.1/2 .1/2 ,1/2 .1/2 900 ^ ' 225 

^0,0 ^0,1 ^1,0 ^1,1 

This implies that in order for the homology computation to be accurate with 
high confidence, we must choose the discretization size M in such a way that 

M'^N'^ for iV ^ 00, 

in contrast to the one-dimensional case. Numerically validated computa- 
tional results which provide more insight into the asymptotic behavior of M 
will be presented in [12]. 

4. Probabilistic tools. The central probabilistic tools used in the re- 
sults of the previous sections are concerned with the asymptotic behav- 
ior of sign-distribution probabilities of parameter-dependent Gaussian ran- 
dom variables. More precisely, let T{5) = (Ti(5), . . . , T„((5))* G R" denote 
a one-parameter family of M"-valued random Gaussian variables over a 
probability space (17,.^, P), indexed by 5 > 0, and choose a sign sequence 
(si, . . . , Sn) G {±1}". We are then interested in the asymptotic behavior of 
the probability 

(20) P{5) = ¥{sjTj{5) > for all j = 1, . . . ,n}. 

To this end, we derive two results. While the first result establishes the 
precise asymptotic behavior of P{S) as 5 — > 0, the second result furnishes an 
upper bound under considerably weakened hypotheses. 
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4.1. Precise asymptotics. We begin with a result which establishes the 
precise asymptotic behavior of sign-distribution probabihties. 

Proposition 4.1. iei (si, . . . G {±1}" denote a fixed sign sequence 
and consider a one-parameter family T{5), 5 > 0, of -valued random 
Gaussian variables over a probability space {Q,J-,¥) which satisfies the fol- 
lowing assumptions: 

(i) For each 6 > 0, assume that the Gaussian random variable T{6) 
has mean G and a positive-definite covariance matrix C{5) G M"^", 
whose positive eigenvalues are given by Xi{5), . . . , A„(5). The corresponding 
orthonormalized eigenvectors are denoted by vi{6), . . . ,Vn{S) . 

(ii) There exists a vector vi = {vn, . . . jVinY G M" such that vi{6) vi 
as 5 ^ and Sj ■ vij > for all j = 1, . . . ,n. 

(iii) The quotient Xi{6)/ Xk{S) converges to as (5 — > for all k = 2, . . . ,n. 

Then the probability P{5) defined in (20) satisfies 



(21) limP((5) 



/detC(5) _ r(n/2) 



<5^o ^ ^ Y Ai((5)" 2 • W2 • (n- 1)! 
Recall that T{l/2) = ■k'^I'^, T{1) = 1 and T{t + 1) = tT{t) for t > 0. 

Proof. Define the diagonal matrix S = {si6ij)ij=i^,,,^n, where 5ij de- 
notes the Kronecker delta, and let = {z G : zj > for j = 1, . . . ,n}. 
Finally, let 

D{6) = Xi{d)-SC{6y^S. 

Using the density of the Gaussian distribution of T{6) according to [3], 
Theorem 30.4, which exists since C{6) is positive-definite, in combination 
with a simple rescaling, the probability in (20) can be rewritten as 



^ ^ VdetC(5) Jz+ \Tr^detC{d) Jz+ 

According to our assumptions, the eigenvalues /xi((5), . . . , /z„(5) of the ma- 
trix D{6) are given by 

Xi (5) 

fii{6) = l and /ifc(^) = for A; = 2, . . . , n, 

Xk{o) 

with corresponding orthonormalized eigenvectors Wk{S) = Svk{5) for k = 
1, . . . , n. Now, let B(6) denote the orthogonal matrix with columns wi{5), . . . , 
Wn{S) and introduce the change of variables z = B{5)C- Moreover, let 



z{Ci,s) = (C2, . . . , Cn) : E CkWkiS) G z+ c : 



on— 1 



k=l 
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and define 

I{CuS)= I exp(-X^//fc(5)C|)d(C2,...,Cn). 

Due to (ii) and the definition of tlie signs Sk, the eigenvector 'Wi{5) has 
strictly positive components for all sufficiently small 5 > and therefore the 
identity 



k=l 



implies that 



(22) 



Jz+ Jb(s)-^z+ \ t^i J 



Z+ JB{S)-^Z+ 

From the definition of I{(^i,6), we can easily deduce that 

i{Ci,s) = cr'- I eM-Cl■jz^'k{5)r|i\d{m 

Jzil,5) \ J 

and this representation implies that 

(23) /(Ci,5)<Cr"^-voln-i(^(l,5)) for ah Ci > and (5 > 0. 

Again, according to (ii), the (n — 1) -dimensional volume of the simplex Z{1, 5) 
converges to the (n — 1) -dimensional volume of the simplex 

Z = {ze Z+:{z- Svi, Svi) = 0} C M", 

which can be computed as 

1 



vol„_i(Z) 



(n-1)! 



Now, let > be arbitrary, but fixed. Note that since we did not make any 
assumptions about the asymptotic behavior of the eigenvectors W2{5),. . . , WniS) 
for (5 — > 0, the sets Z{1,6) do not have to converge. Yet, (ii) yields the exis- 
tence of a compact subset K C R""^ such that Z(l,6) <Z K for all sufficiently 
small 6 > 0. Due to (iii), the integrand in (22) converges to 1 uniformly on K 
and we therefore have 

lim /(Ci , 5) = Ci~^ ■ vol„„i (Z) for all Ci > 0. 

5— >0 
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Due to (23) and vol„_i(Z(l, 5)) ^ vol.„_i(Z), we can now apply the domi- 
nated convergence theorem to pass to the hmit (5 — > in (22), which leads 
to 

hm / e-^*^(^)^ dz = voL„_i(Z) • Cr'^e-^? dCi = vol„_i(Z) • 

&^^Jz+ Jo 2 

This completes the proof of the proposition. □ 

The above result is the main probabilistic tool for our one-dimensional 
results in Section 2. In this situation, Proposition 4.1 is used to study random 
variables in M^. We note that for three-dimensional random variables, it is, 
in fact, possible to explicitly compute the probability P{6) for each 6 > 
as a function of the entries in the covariance matrix C{5). This has been 
demonstrated in [15] using the inversion formula ([5], Formula (29.3)) and 
could be used to study the asymptotic behavior of P{5) as 5^0. However, 
for the results of Section 3, we must also consider higher-dimensional random 
vectors and the method in [15] no longer applies. 

4.2. Upper bounds. While the above result establishes the precise asymp- 
totic behavior of the sign-distribution probability P{S), the necessary as- 
sumptions for Proposition 4.1 are not satisfied in all instances discussed in 
Section 3. Yet, in these cases, it suffices to obtain an upper bound on the 
probability P{S). Such an upper bound is the subject of the following result, 
which can be derived under fairly weak assumptions. 

Proposition 4.2. Let (si, . . . , Sn) G {±1}"' denote a fixed sign sequence 
and consider a one-parameter family T{5), S > 0, of -valued random 
Gaussian variables over a probability space (r2,JF, P) which satisfies the fol- 
lowing assumptions: 

(i) For each 6 > 0, the Gaussian random variable T{6) has vanishing 
mean G M" and a positive-definite covariance matrix C{5) G M"^" whose 
positive eigenvalues are given by Ai(5), . . . , A„((5). Corresponding orthonor- 
malized eigenvectors are denoted by vi{6), . . . ,f„((5). 

(ii) There exists a vector vi = {vn, . . . ,vinY G such that vi{5) vi 
as 5 — > and sj ■ vij > for all j = 1, . . . , n. 

Then the probability P{5) defined in (20) satisfies 



n 
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Proof. We adopt the notation used in the proof of Proposition 4.1. We 
then have 



is a consequence of (ii), and this immediately imphes (24). □ 

Note that Proposition 4.2 does not require any assumption of the eigen- 
value Ai(5). In particular, Xi{6) need not be the smallest positive eigenvalue 
of C{6), in contrast to Proposition 4.1. We also note that we could easily 
modify the above result to yield lower bounds for the asymptotic behavior 
of P{S). Since these are not needed for our applications, we refrain from 
formulating them here. 
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